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I.  INTRODUCTION 

Under  a large  positive  axial  acceleration  field,  a rocket  will 
compress  along  its  length  similar  to  a shaft  under  Eulerian  loading.  The 
strain  energy  stored  will  cause  a return  to  the  undeformed  configuration 
if  the  acceleration  decay  is  sufficiently  smooth  (Figure  1).  If,  however, 
the  acceleration  is  rapid,  as  shown  in  Figure  2,  an  oscillation  will 
result  about  the  abcissa.  In  these  figures  the  acceleration  profiles 
are  shown  in  the  insert  while  the  corresponding  displacements  appear  in 
the  main  figures.  If  a rocket  with  a steep  acceleration  decay  also 
possesses  a flexible  or  loose  joint,  it  is  conceivable  that  its  snap-back 
response  can  be  destabilizing. 

The  VIPER  has  all  of  the  characteristics  previously  cited.  It  has 
a very  high  peak  acceleration,  a fairly  steep  acceleration  decay,  and  a 
joint  which  appears  to  have  little  stiffness  in  the  axial  downrange 
direction.  Additionally,  flight  tests  have  indicated  a relatively  high 
rocket-to-rocket  dispersion.  All  of  these  characteristics  led  the  VIPER 
Project  Office  to  request  a snap-back  analysis.  This  report  describes 
the  approach  used  to  model  this  event  and  presents  some  of  the  results. 


II.  APPROACH 

The  analysis  of  the  VIPER  snap-back  phenomenon  was  conducted 
in  two  phases.  The  purpose  of  these  analyses  was  to  determine  the 
snap-back  forces  in  the  magneform  joint.  Tests  performed  earlier  indi- 
cated that  forces  of  5000  lb  or  greater  would  cause  joint  separation 
which,  in  turn,  may  result  in  flight  destabilization.  Because  the  stiff- 
ness of  the  joint  in  tension  could  not  readily  be  calculated,  it  was 
assumed  to  be  the  same  as  in  compression.  Then,  knowing  the  displacement 
at  the  joint,  the  force  imposed  on  the  joint  could  be  computed.  First, 
a detailed  finite  element  static  analysis  was  made  to  determine  deflection 
versus  length  produced  by  the  launch  acceleration  of  7000  g imposed  as 
a static  load.  This  model  is  presented  in  Figure  3.  Next,  a lumped 
parameter  model  was  formulated  for  the  dynamic  analysis  (Figure  4)  and 
checked  against  the  more  complex  static  model  by  imposing  the  launch 
acceleration  statically  and  comparing  the  resulting  deflection  field 
with  that  derived  using  the  finite  element  procedure.  Finally,  the 
acceleration  curve  was  approximated  by  a piecewise  linear  function  and 
was  imposed  on  the  dynamic  model. 


III.  STATIC  ANALYSIS 


The  static  analysis  employed  a finite  element  program,  AMG054, 
which  analyzes  stresses  and  deflections  due  to  symmetric  loading  on  a 
body  of  revolution.  The  program  is  capable  of  accommodating  a large 
number  of  degrees  of  freedom  (in  this  case,  2160).  AMG054  predicts  the 
stress  and  displacement  due  to  bending  and  axial  loading.  Figures  3,  5, 
6 and  7 show  the  AMG054  model  and  the  warhead  details.  The  loading 
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was  defined  as  a distributed  axial  load  with  a magnitude 


7000  g I m(x)dx 


and  a fore-to-aft  direction.  Figure  8 shows  the  deformed  configuration 
resulting  from  the  7000-X  386-lb  force.  Displacements  were  multiplied 
by  100  for  "visualization." 


IV.  DYNAMIC  ANALYSIS 

The  dynamic  model  (Figure  4)  is  a lumped  parameter  system 
where  the  mass  of  the  structure  is  concentrated  at  discrete  points 
along  the  length,  and  its  stiffness  is  represented  by  massless  rods. 

The  equations  of  motion  of  this  structure  form  a coupled  system  of  five 
second-order  differential  equations  with  constant  coefficients.  This 
system  may  be  written  symbolically  as 

mijXj  + kij  Xm  " Fi(t)  i,j  = 1,  2 5 . (2) 


The  stiffness  of  the  ith  rod  is  defined  as 


*i  Lt  » 

where  is  the  cross-sectional  area  of  the  rod,  is  its  length,  and 
is  the  Young's  modulus  of  the  material. 

Next,  the  equations  of  motion  are  written  as  the  system  of  10  first- 
order  differential  equations  as  follows: 

DV(2)/Dt  - V(3) 


DV(3)Dt  - ^-["k1V(2)  + kjV (4)J  + g(t) 


DV(4)/Dt  - V(5) 


DV(5)/Dt  - MV(2)  ’ <kl  + k2>V<3)  “ k2V(6)]  + g(t) 
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DV(6)/Dt  = V(7) 

DV(7)/Dt  = ^~[k2V(4)  - (k2  + k3)V(6)  + k3V(8)]  + g(t) 

DV(8)/Dt  = V(9) 

DV(9)/Dt  - [k3V(6)  - (k3  + k4)V(8)  + k4V(9)]  + g(t) 

DV(10)/Dt  - V(ll) 

DV(ll)/Dt  = ^[k4V(8)  - (k4  + k5)V(10)]  + g(t)  , (4) 

where  g(t)  is  the  acceleration  profile  described  in  Figure  9,  and 
V(l)  = t = Time 


V(2)  - Xj_ 

V(3)  = xL 
V(4)  = x2 
V(5)  = x2 

V(6)  = x3 

V(7)  = ^ 

V(8)  = x4 
V(9)  - x4 
V(10)  • x5 

V(ll)  - x5  . (5) 
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Computer  programs  INTERP  and  SSIMDE  were  used.  INTERP  performs  an 
interpolation  of  the  piecewise  linear  functions  g(t)  for  each  time 
step  in  the  integration  procedure.  SSIMDE  employs  the  Runga  Kutta 
technique  to  integrate  Equation  (4) . The  program  used  for  data  input 
and  graphics  is  shown  in  the  appendix. 

To  check  this  procedure,  a five-lumped  parameter  model  of  a uniform 
beam  was  constructed.  Figure  10  shows  the  model  and  significant  para- 
parameters,  while  Figures  11  through  15  show  the  displacements  for 
m^  through  m^.  The  overall  displacement  computed  was  0.0454  in., 

while  the  closed  form  solution  was  0.0438  in.  Because  the  error  was 
only  3.67„,  the  model  appears  accurate. 

The  VIPER  warhead  model  masses  and  stiffnesses  were  calculated  and 
a trial  run  was  made.  A run  was  made  with  g(t)  increasing  to 

2.702  X 10^  lb  and  remaining  at  that  level  until  the  displacements  of 
the  model  reached  a stable  value  for  all  masses.  The  resulting  displace- 
ments were  compared  to  the  finite  element  predictions  and  corrected 
accordingly. 

Finally,  four  acceleration  profiles  were  considered; 

a)  Case  1 - g(t)  decays  linearly  from  0.0055  to  0.03  sec. 

b)  Case  2 - g(t)  decays  linearly  from  0.0055  to  0.01  sec. 

c)  Case  3 - g(t)  decays  linearly  from  0.0055  to  0.006  sec. 

d)  Casre  4 - g(t)  decays  linearly  from  0.0055  to  an  acceleration 

of  1250  g at  0.008  sec  and  then  to  0 at  0.01  sec. 


V.  RESULTS  AND  CONCLUSIONS 

The  displacements  resulting  from  the  four  cases  are  shown  in 
Figures  16  through  19.  Cases  1,  2,  and  3 were  performed  to  investigate 
the  sensitivity  of  the  snap-back  forces  and  displacements  to  the  slope 
of  the  decay  portion  of  the  acceleration  curve.  Case  4 is  a fairly 
accurate  representation  of  the  actual  acceleration  profile.  Table  1 
presents  the  maximum  negative  displacements  at  mass  number  five  which 
occur  at  snap-back  in  each  case  and  the  corresponding  force  in  the  joint. 


Sellers,  W.  R. , Jr.  and  Gibbs,  B.  G.,  Descr ip tions-General- Purpose 
Computer  Subroutines.  US  Army  Missile  Command,  Redstone  Arsenal,  Alabama, 
January  1975,  Addendum  - January  1977,  Report  No.  TR-WS-75-2. 
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TABLE  1.  MAXIMUM  NEGATIVE  DISPLACEMENTS 
AT  MASS  NUMBER  5. 


Case 

-X^ (Max) 

(Max) 

1 

-0.00007 

800.0 

2 

-0.00076 

987.5 

3 

-0.00142 

1851.0 

4 

-0.00101 

1322.1 

The  obvious  conclusion  is  that  the  snap-back  force  increases  as  the  slope 
of  the  acceleration  decay  curve  increases.  Another  run  was  made  at  one- 
half  the  period  of  the  fundamental  longitudinal  natural  frequency.  This 
resulted  in  an  extremely  large  deflection  (-0.01010  in.)  and  force  on  the 
joint  (13,130  lb). 

Case  4 is  a realistic  representation  of  the  actual  acceleration 
profile.  The  predicted  load  at  the  joint  is  lower  than  5000  lb  and  is, 
therefore, of  insufficient  magnitude  to  cause  separation.  However,  it 
must  be  cautioned  that  if  the  decay  period  approaches  one  half  of  the 
fundamental  longitudinal  frequency,  the  excursion  and,  consequently, 
the  force,  will  become  very  high. 

Still  another  factor  which  should  be  considered  in  future  work 
is  the  motor  case  weight.  Since  the  motor  case  is  approximately  equal 
to  the  forward  portion  of  the  rocket  a model  of  the  entire  rocket  should 
be  developed  and  tested.  Unfortunately  there  was  not  time  to  accomplish 
this  during  this  effort. 
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Figure  3.  Finite  element  model  (static) 


Figure  5.  Detail  of  aft  warhead  body 
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Figure  6.  Detail  of  warhead 


Figure  8.  Deformed  warhead 
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Figure  9.  Rocket  launch  acceleration. 
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Appendix.  INPUT  PROGRAM 


PROGRAM  MAIN 


74/74 


OPT»l 


FTN  4.2*74355 


1 0 


lc 


20 


?5 


39 


35 


40 


45 


50 


55 


PROGRAM  MA 1 N II NPUT *A5 , OUTPUT*65  . TAPE6=0UTP(JT  « 7 APF2=5 1 3 . 
1TAPF5-65) 

DIMENSION  VT(lt) • W ( 3 > 

DIMENSION  X(1000) .Y(IOOO) . TXST (4) , IYST ( 14) ,INUM(5> 
COMMON/FORC/XT (6) ,YT(6) 

FXTFRNAL  DESUB 

0A7A  lXST/84.73.77,69/ 

DATA  I YST/68. 73, 83. 80 *76. 65 .67, 69, 77, 49,78. 84, 3?, 3?/ 
DATA  INUM/4Q.50.51 .52.53/ 

DATA  XT/O. 0. 0.0025. 0.0055. 0.008.. 81.1,/ 

DATA  YT/O. 0.6500.. 7000.. 1250. ,0.»0./ 

REMIND  2 
TMAXxO.O? 

NP»1 

J*0 

IC=0 

DEL*. 00001 
DELmN*.5*0E| 

WRITE (6.1000)  (I, I. 1=1. 5) 

DO  100  J=l.ll 
100  VT ( I ) =0 .0 

WRITE (2)  VT(1)»VT(2)»VT(4>»VT(6)»VT(8)»VT(10> 

WRITE (6.1 1 09)  VT(1) . (VTII) .1*2,10.2)  . 

1 (VT (K) ,K*3. 1 1 .2) 

200  J*J»1 

CALL  SSIMDE (VT.W.DEL .DFLMN. 10. IC.DF5UR) 

IF (MOD ( J. 1 0 ) .NE.  0)  60  TO  150 
WRITE (6.1 100)  VT(1).(VT (I) ,1*2.10,2). 

1 ( VT (K ) .K=3 , 1 1 .2) 

NP=NP» 1 

WRITE (2)  VT (1) .VT (2) ,VT(4) .VT(6) .VT (8) .VT (10) 

150  CONTINIJF 

IF « VT  <11  ,GT.  TMAX)  60  TO  250 
GO  TO  200 
250  CALL  INITTM) 

DO  400  K=l«5 
IYST(14)=INHM(K) 

CALL  RINITT 
REWIND  2 
DO  300  m*1,mP 

READ (2)  X(M)»VT(1),VT(2).VT(3),VT(4)»VT(5) 

300  Y(M)*VT(K) 

CALL  NPTS(NP) 

CALL  XERM (4 ) 

CALL  TERM (4) 

CALL  CHECK (x.Y) 

CALL  DSPLAY(X.Y) 

CALL  M0VA8S (20.550) 

CALL  VLARFL (14.IY5T) 

CALL  NOTATE (500.20.4, IXST) 

CALL  movEA(o. 0.0.0) 

CALL  DRAWA(TMAX.O.O) 

CALL  SLIMX (=50,900) 

CALL  SL IMY ( 450,700 ) 

CALL  NPT5 (6) 

CALL  DLIMxn.O.TMAX) 

CALL  OL!MY(p. 0.10000.) 
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PROGRAM  MAIN 


74/74  OPT=l 


FTN  4.2*74,355 


CALL  CHECK (XT. YT> 

CALL  05PLAY (XT »YT) 

CALL  SLJMX ( 1 50*400) 

CALL  SL IMY (125,700) 

CALL  BELL 
CALL  TINPUT ( ICR) 

CALL  ERASE 
400  CONTINUE 

1000  FORMAT ( 1H1 * TIME  »5(*  X(«I1*)/X  OOT(*»1*)  •)) 
1100  FORMAT ( 1H  F«.5*5F15.7/QX,5F15.7) 

eno 


SUBROUTINE  OESIIB 


74/7* 


OPT  = 1 


FTN  *.2*7*355 


SUBPOUT 1 ME  nESUR(VT.F.J) 

COMMON/FOPC/XT (6) « VT  (6) 

OIMFNSION  VT ( 1 ) 

CONS*l.EB 

5 0*0.0 

IF (MODI J»2>  .FO.  0)  CALL  INTFPP < VT ( 1 ) , XT.TT»6,2.G.NFPP> 
6*0*386.* 

GO  T0<10. 20. 30. *0.50. 60. 70. 80. 90. 100) . J 
10  F*vT (3) 

10  GO  TO  200 

20  F=-(6*.00*<VT  <2>-VT<*> ) >*CONS*G 
GO  TO  200 
30  F=VT (5) 

GO  TO  200 

15  *0  F=( 32 .*VT(2) -80.00* VT<4) *47.00*VT (ft) ) *CONS*G 

GO  TO  200 
50  F*VT  < 7 ) 

GO  TO  200 

60  F*(13.00*VT (*) -2*.0*VT (6)*11.*VT(8)) *CONS*G 
20  GO  TO  200 

70  F = VT ( 9 ) 

GO  TO  200 

80  F*  (7.00*VT  (M-19.00*VT  (8)  *13.0*VT  ( 10)  ) *rONS*G 
GO  TO  200 

25  oi)  F = VT  111) 

GO  TO  200 

100  F=(47.00*VT(8)-67.000*VT(10) )*CONS*G 
C 200  WPTTE(6.100O)J,F.G 
200  RETIIPN 

■»0  1000  FOPMATC1H  * J s*I5.5X*  F * *E15.7.5X*  fi  = *E15.7) 

END 
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